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PREFACE 
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THERMODYNAMICS  OF  FRICTIONAL  MATERIALS 


A  CONSTITUTIVE  THEORY  OF  SOILS  WITH  DILATANT  CAPABILITY 

PART  I:  INTRODUCTION 

1.  The  main  task  of  this  work  is  to  use  the  concepts  of  endochronic  plasticity  to  de¬ 
rive,  from  thermodynamic  principles,  a  constitutive  equation  for  soils  with  dilatant  capability. 
The  aim  has  been  to  set  the  theory  in  a  three  dimensional  framework  and  endow  the  resulting 
constitutive  relations  with  the  capability  for  accounting  for  the  effects  of  complex  deformation 
histories.  The  result  is  a  domain  of  applicability  broad  enough  to  capture  sufficiently  complex 
phenomena  to  make  the  constitutive  equations  useful  in  applications  to  realistic  deformation 
histories  encountered  in  practice. 

Overview  of  Endochronic  Plasticity 

2.  The  endochronic  theory  (Valanis,  1971;  Valanis,  1980;  Valanis  and  Read,  1980;  Valanis 
and  Lee,  1982;  and  Valanis,  1986)  deals  with  the  plastic  response  of  materials  by  means  of  memory 
integrals  expressed  in  terms  of  memory  kernels  which  contain  the  essence  of  the  material  behavior. 
The  kernels  are  invariably  fast  decaying  functions  indicating  a  material  memory  which  fades  at 
a  high  rate.  In  the  specific  case  where  these  kernels  are  approximated  by  Dirac  delta  functions 
(which  decay  at  an  infinite  rate)  one  recovers  a  theory  with  yield  surfaces.  Thus  yield  surface 
plasticity  is  a  subset  of  the  endochronic  theory. 

3.  The  material  memory  is  specified  with  respect  to  an  intrinsic  time  which  is  the  path 
traversed  by  the  plastic  strain  in  a  plastic  strain  space.  If  the  material  is  isotropic  in  its  reference 
configuration  the  path  length  is  given  in  terms  of  one  material  parameter.  See  Ekiuation  (5),  Part 

n. 

4.  The  application  of  endochronic  plasticity  to  soils,  insofar  as  work  that  has  appeared 
in  the  literature  is  concerned,  has  been  limited  either  to  one  dimensional  histories  (shear  and 
hydrostatic  response)  under  cyclic  and  other  loading-unloading  conditions  or  to  two  dimensional 
histories  (shear-volume  interaction)  under  densifying  conditions.  The  question  of  dilatant  behav¬ 
ior  had  hitherto  remained  unresolved.  The  reasoning  that  led  to  the  present  treatment  is  broadly 
as  follows; 

The  coupling  between  deviatoric  and  hydrostatic  behavior  that  will  ultimately  lead  to 
dilatant  deformation  must  come  from  three  sources: 
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a.  The  intrinsic  time  through  the  parameter  k  in  Ekiuation  (5); 

b.  The  expression  for  the  free  energy; 

c.  The  rate  equations  for  the  internal  variables. 

5.  Source  a  alone,  will  always  give  rise  to  densification  as  the  application  of  the  relevant 
equations  to  concrete  by  Valanis  and  Read  (1986)  actually  showed.  Source  b  is  not  physically 
viable  because  given  a  soil  with  a  certain  porosity  the  onset  of  dilatancy  under  monotonic  shearing 
is  governed  by  the  prevailing  hydrostatic  stress.  If  b  is  to  be  the  source,  then  a  change  in  the  form 
of  the  free  energy  must  take  place  upon  varying  the  hydrostatic  stress  which  is  not  physically 
plausible.  In  other  words,  one  would  expect  the  form  of  the  free  energy  to  remain  invariant  with 
a  change  in  the  hydrostatic  stress.  The  remaining  plausible  cause  is  source  c  and  this  is  the  one 
that  we  developed  in  the  present  work. 

Overview  of  Report 

6.  In  Part  II,  a  simple  model  is  developed  to  illustrate  basic  principles.  The  point  of 
departure  is  the  critical  state  theory  (SchoBeld  and  Wroth,  1968)  which  gives  physically  valid 
predictions  in  stress  states  near  equilibrium.  The  ideas  of  critical  state  are  developed  beginning 
on  page  6  where  a  simple  constitutive  equation  is  given  which  exhibits  the  salient  features  of  den- 
sifying  and  dilatant  behavior.  In  particular,  the  groundwork  is  laid  with  regard  to  the  physical 
properties  of  the  coupling  coefficients  that  allow  for  prediction  of  dilatancy.  The  simple  constitu¬ 
tive  equation  developed  Part  II  is  limited  in  a  number  of  ways,  not  the  least  of  which  is  the  fact 
that  it  is  constrained  to  obey  an  isotropic  hardening  law. 

7.  An  isotropic  hardening  law  may  do  a  reasonable  job  if  no  unloadings  or  sharp  kinks 
are  encountered  in  the  deformation  path,  but  otherwise  it  will  give  poor  results.  What  does  not 
seem  to  have  been  realized  in  the  community  of  practitioners  and  other  plasticians  is  that  the 
essence  of  plastic  behavior  in  materials  (almost  without  exception)  insofar  as  the  shear  response  is 
concerned  is  kinematic  hardening.  In  practical  thermodynamic  terms  this  means  more  than  one 
internal  variable  (in  fact  a  large  number  or  spectrum)  is  needed  if  one  is  to  obtain  a  realistic  shear 
response.  This  possibility  is  explored  in  Part  III  where  the  foundations  are  laid  for  two  distinct 
approaches,  Theory  I  on  page  20  and  Theory  II  on  page  28,  where  two  different  hydrostatic 
constitutive  equations  are  derived.  Theory  II  is  found  to  be  preferable  in  the  sense  that  it  is  not 
subject  to  analytical  difficulties  and  approximations  and  is  physically  more  appealing.  In  Part 
IV  a  specific  application  is  presented  using  a  special  form  of  Theory  II. 

8.  In  conclusion,  the  analytical,  physical  and  thermodynamic  foundations  are  laid  for 
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the  derivation  of  constitutive  equations  that  can  account  for  densifying  or  dilatant  behavior  of 
soils  depending  on  the  prevailing  density,  hydrostatic  stress  and  deformation  history.  Specifically, 
it  is  believed  that  Equations  (152)  and  (153)  have  the  proper  foundation  and  appropriate  ana¬ 
lytical  form  to  deal  with  complex  forms  of  soil  behavior.  The  application  of  these  equations  to 
deformation  histories  encountered  in  practice  will  be  the  subject  of  further  research. 
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PART  II:  FRAMEWORK  OF  A  SIMPLE  MODEL 


9.  This  work  begins  with  the  simplest  model  of  coupled  linear  relations  between  stress 
and  strain-rates.  The  rates  are  with  respect  to  the  intrinsic  time  scale  z  and  the  coupling  is 
between  the  hydrostatic  and  deviatoric  strain  rates  on  one  hand  and  corresponding  stress  states 
on  the  other.  Specifically,  we  begin  with  the  Equations  (1)  and  (2)  as  follows: 


de'* 


dfP 


8  =  aii-j-  +  312-3— 
dz  dz 


(1) 


a  =  321 


deP  dfP 


(2) 


in  the  usual  notation  where  bold  symbols  denote  tensor  quantities  and  a  dot  between  two  tensors 
denotes  their  scalar  product.  The  superscript  p  implies  that  the  strain  rates  are  pleustic  as  defined 
in  terms  of  the  total  shear  and  hydrostatic  components  as 


deP 


de  -  -^ds 

2p 


dfP  =  de  -  —da 
K 


(3) 

(4) 


where  /i  and  K  are  the  elaistic  shear  and  bulk  moduli  respectively.  The  intrinsic  time  meeisure  dz 
is  given  by  Equation  (5) 


dz*  =  lldePll*-»-fc*ld€Pl*  (5) 

where  ||  •  j|  denotes  the  norm  of  the  tensor  and  k  is  the  coupling  parameter  between  the  deviatoric 
and  hydrostatic  deformation.  Equations  (1)  and  (2)  are  the  classical  thermodynamic  rate  equa¬ 
tions  where  the  “forces”  s  and  a  on  one  hand  and  the  “fluxes”  deP  jdz  and  dt^  jdz  on  the  other 
are  related  by  the  matrix  of  “Onsager”  coefficients 

On  ai2 
321  O22 

which  traditionally  have  obeyed  the  symmetry  condition 


(6) 


3r»  —  3*r 


(7) 
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However,  as  will  be  discussed  presently,  the  symmetry  condition  is  not  essential  so  long  as  the 
positive  dissipation  inequality  (8)  is  satisfied.  Thus,  no  prior  assumptions  are  made  about  the 
symmetry  of  ar*  with  respect  to  r  and  a  and  the  possibility  that  au  7^  821  wil!  be  considered. 
Note  that  both  of  these  are  deviatoric  because  of  Equations  (1)  and  (2).  In  effect  ai2  is  deviatoric 
because  s  is  likewise,  and  only  the  deviatoric  part  of  321  plays  a  role  in  the  inner  product  on  the 
right  hand  side  of  Equation  (2).  However,  the  positive  rate  of  dissipation  must  be  observed  in 
which  case: 


deP  dcP 

s-  +  <7—  >  0 
dz  dz 


(8) 


for  ((dePjl  >  0  and/or  \de'’\  >  0.  It  follows  from  inequality  (8)  and  Equations  (1)  and  (2)  that 


an 


dz 


,  ,  dePdfP  (deP\^  ^ 

+  +  +  j  >0 


10.  Consider  the  constraints  on  the  matrix  [A]  where 


(9) 


flu  ai2 
321  022 


(10) 


which  result  from  inequality  (9).  To  this  end  set 


312  +  321  =  2B 


(11) 


and  let 


B- 


deP 

dz 


dfiP 

— r-p  COB 
dz 


(12) 


where  =  ||B||  and  deP  =  HdePl]. 

11.  It  follows  that  inequality  (9)  will  not  be  violated  if  the  discriminant  of  the  quadratic 
form  Q  given  by  Equation  (13). 


Q  =  an 


deP 


dz 


+  2)3  cos  ^ 


deP 

dz 


dcP 


dz 


+  022 


(f)' 


(13) 


is  negative  and  on  and  022  are  positive.  Thus  the  necessary  and  sufficient  condition  is 


an  >  0,  022  >  0»  /3*cos*^—  011022  <  0  (14) 

for  all  <j>.  The  worst  possible  czise  is  when  ^  =  0.  Thus  the  necessary  and  sufficient  condition  that 
inequality  (9)  be  observed  is  that 
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Figure  1:  Hydrostatic  response  depicted  by  Equation  (18). 


an  >  0,  022  >  0,  aiia22  -  |  ||ai2  +  a2il|*  >  0  (15) 

since 

^*  =  ||Bl|^  =  l|l(at2  +  a2,)||*  (16) 

Physical  Identification  of  Matrix  [A] 

12.  In  this  section  connections  are  established  between  the  proposed  simple  constitutive 
theory  specihed  by  Ekiuations  (l)  and  (2),  and  other  existing  theories  which  have  proved  physically 
sound  in  describing  qualitatively  the  constitutive  response  of  soils  under  monotonic  loading.  One 
such  theory  is  that  of  the  critical  state  developed  by  Roscoe  et.  al.  (1958).  The  fundamental 
assumption  that  underlies  this  theory  is  that  the  rate  of  plastic  work  is  proportional  to  the 
hydrostatic  stress.  This  statement  is  given  analytically  in  incremental  form  by  Equation  (17). 

8  •  de*’ +  crdc*’ =  aMde*’  (17) 

One  observes  that  under  purely  hydrostatic  conditions  Equation  (17)  leads  to  the  result 

adcP  =  0  (18) 
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which  implies  that  either  =  0  or  =  0,  that  is,  either  a  =  0  or  the  response  is  elastic. 
Equation  (18),  therefore,  leads  to  the  hydrostatic  response  depicted  in  Figure  1  when  the  elastic 
part  of  the  response  is  linear.  Also  depicted  in  that  figure  is  the  actual  response  which  qualitatively 
looks  like  the  critical  state  model.  In  presentation  of  the  criticzd  state  theory  (e.g.  Schofield  and 
Wroth,  1968),  liquation  (17)  is  usually  restricted  to  the  case  where  de**  ^  0.  However,  one  would 
rather  have  the  response  for  all  loading  conditions  to  be  contained  within  the  theory  without 
additional  stipulations.  ^  Since,  otherwise,  the  critical  state  theory  gives  physically  sound  results, 
Elquations  (1)  and  (2)  must  contain  Equation  (17)  in  some  sense.  Specifically  Equation  (2)  could, 
in  fact,  contain  Elquation  (17)  if 

aji  =  fljis,  022  =  (19) 

where  and  022  are  constants.  In  this  specific  instance  Equation  (28)  would  have  the  following 
incremental  form: 

adz  =  aJiS  •  de**  +  o52ffdc**  (20) 

In  the  very  particular  case  where  a^x  =  <>22  Equation  (20)  would  read 

adz  =  022(9  •  de^  -f  ode**)  (21) 

which  is  a  statement  to  the  fact  that  the  increment  in  plastic  work  is  proportional  to  the  hydro¬ 
static  stress,  except  that  dz  is  not  equal  to  dt^  as  it  would  be  in  the  case  of  the  critical  state 
theory.  However,  one  would  like  to  retain  the  greater  freedom  of  Equation  (20)  rather  than  bound 
by  the  Equation  (21).  Again  observe  that  under  purely  hydrostatic  conditions  and  in  view  of 
E^quations  (8)  and  (20)  one  finds  that 

ad€^{k  -  o$2)  =  0  (22) 

Now  if  k  ^  022  again  obtains  the  hydrostatic  response  in  accordance  with  the  critical  state 
theory,  given  by  Equation  (18)  and  depicted  in  Figure  1.  However,  there  is  another  choice.  By 
setting 

'The  indeterminate  nature  of  the  response  for  purely  hydrostatic  loading  was  the  motivation  of  introducing  the 
elliptical  surface  into  the  critical  state  theory  by  Burland  (196S)  through  the  introduction  of  an  alternative  form 
of  the  right-hand  side  to  Equation  (17).  Equation  (17),  which  is  based  on  the  Taylor  energy  relationship,  typically 
fits  experimental  data  well  and  rightly  held  an  axiomatic  position  in  the  original  theory.  Here,  the  necessity  for 
stipulating  de'  0  will  be  removed.  The  more  general  problem  of  indeterminate  response  under  hydrostatic  loading 
wiU  be  addressed  beginning  on  page  IS. 


Equation  (22)  is  satisfied  identically  so  that  nothing  can  be  said  about  either  a  or  and  any 
purely  hydrostatic  response  is  admissible  (such  as  the  one  found  in  experiment). 

13.  In  view  of  the  above  discussion  set: 


an  —  OnS 
aji  = 

ajj  = 


in  which  case  Equations  (1)  and  (2)  become 


(24) 

(25) 

(26) 


(27) 


a 


deP 


=  0218.  — +  a„o  — 


(28) 


Constitutive  Response  for  Case  of  0J2  =0 

14.  This  section  examines  the  role  of  the  coefficient  o^  in  Equation  (27).  It  follows  from 
this  equation  that 


Oil 


d^ 


dz 


dcP 


(29) 


Whatever  the  sign  of  of]  a  case  can  be  made  that  the  shear  response  depicted  by  Equation  (29) 
is  unrealistic.  For  instance  in  the  light  of  the  convention  that  is  positive  in  compression  let  0^2 
be  positive.  One  can  then  envision  a  situation,  where  the  hydrostatic  deformation  is  densifying, 
in  which  case  the  shear  response  will  become  unstable  if 


“12 


dcP 


dz 


Since  such  possibilities  are  not  physically  realistic  set  a?.  =  0. 


(30) 


Existence  of  a  Yield  Surface 

15.  In  this  section  Equations  (27)  and  (28)  are  used  in  conjunction  with  Equation  (23) 
and  the  definition  of  the  intrinsic  time  scale,  given  by  Ek)uations  (5)  and  (8),  to  derive  a  yield 
surface.  Substitution  of  the  former  into  the  latter  leads  to  the  following  result: 

—  2dii<22j<r  +  ||o|is||*  =  0  (31) 

The  geometry  of  the  yield  surface  will  depend  on  the  nature  of  the  functional  dependence  of  aix 
on  cr  and  z.  To  examine  various  forms  of  its  geometry  let  ||s||  be  denoted  by  q  in  which  event 
Equation  (31)  becomes 

<7^  -  2aiia2i<r  +  (ofi)*?*  =  0  (32) 

16.  First  examine  the  specific  case  where 

dll  =  aii<Tf(z)  (33) 

In  view  of  Equation  (27)  and  the  fact  that  0^2  =  0)  Equation  (33)  implies  that  the  shear  stress  is 
directly  proportional  to  the  hydrostatic  pressure,  in  accordance  with  the  classical  law  of  friction. 
In  this  case.  Equation  (32)  becomes 

-  1)  =  (34) 

or 


Thus  the  yield  locus  becomes  a  straight  line,  the  slope  of  which  increases  as  the  material  hardens, 
i.  e.  as  f{z)  increases  with  deformation.  When  /(z)  becomes  a  constant  then  a  “failure  line” 
(surface)  is  reached. 

17.  Another  case  worth  investigation  is  the  one  where  on  is  independent  of  a.  This  is 
not  realistic  in  the  case  of  soils  but  is  of  academic  interest  none-the-less.  In  this  case 

an  =  «n/(z)  (36) 

and  the  yield  equation  given  by  Equation  (32)  now  becomes 
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a^-2a5ia?i/(zy  +  (a|,)V  =  0  (37) 

or 

(T  -  a|, at, /(*))’  +  (os,)"?’  =  (at, oi, /(.))’  (38) 

which  is  the  equation  of  an  ellipse  with  center  (a2i®ii/('®^)>0)-  These  two  simple  cases  are  illus¬ 
trated  in  Figure  2. 

18.  In  reality  the  dependence  of  an  on  a  is  likely  to  lie  between  these  two  extremes. 
Thus  with  reference  to  Equation  (32),  an  is  likely  to  be  of  a  form  such  as: 

an  =  <rA{ff)alif{z)  (39) 

where  A{a)  is  decreasing  function  of  <t,  >1(0)  =  1,  and  is  a  monotonically  increasing  func¬ 

tion  of  z,  becoming  constant  as  z  — »  oo.  To  examine  the  effect  of  Equation  (39)  on  the  geometry 
of  the  yield  surface  one  can  choose,  typically,  A(<t)  to  be  the  following  decaying  exponential, 

A{a)  =  (40) 

in  which  case  Equation  (32)  becomes 

cj\2a\J{z)a\,e-<^  -  1)  =  (a^,)^^*  (41) 
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q 


Figure  3:  Yield  surfaces  for  various  values  of  m. 

Now  observe  that  q  is  zero  at  <r  =  0  but  also  at  a  limiting  hydrostatic  stress  a  =  <f>t  such  that 

=  1  (42) 


or 


a^e  =  logm(-r)  (43) 

where  we  have  set  2aiia2if{z)  =  »w(z).  Thus  in  terms  of  these  parameters,  where  m(z)  is  a 
monotonically  increasing  function  of  z  with  the  property 


lim  m(z)  =  moo 

Equation  (41)  becomes 


(44) 


<lo 


(45) 


where  Oji  =  <^el<lo-  Equation  (45)  gives  a  family  of  yield  surfaces  each  corresponding  to  an 
increasing  value  of  m(z).  In  Figure  3,  a  subset  of  these  surfaces  is  shown  for  various  values  of 
m(z).  The  situation  is  more  complex  when  an  is  also  a  function  of  density. 


11 


Strain  Response  Given  the  Stress  History 


19.  Recall  Equations  (27)  and  (28)  which,  in  light  of  the  stipulation  that  =  0,  become 


8  =  an 


dz 


(46) 


Use  of  Equation  (46)  in  (47)  gives  the  relation 

— )  dz  =  a\'i<jd^  (48) 

an  J 

Equation  (48)  is  central  in  determining  whether  the  volumetric  response  will  be  contractive  or 
dilatemt.  Adopting  the  convention  that  cr  and  £**  are  both  positive  in  compression,  then  de^  >  0 
implies  contraction  and  d^  <  0  implies  dilantancy.  Thus  if 


£|i^2 

an 


>  a 


(49) 


then  the  current  deformation  is  dilatant.  FVom  Equation  (32)  it  is  apparent  that  a  and  q  are 
related  since  they  both  have  to  lie  on  the  yield  surface.  Specifically,  if  Equation  (39)  is  substituted 
into  Equation  (32)  the  yield  relationship  is  given  by  the  following: 

a*  -  2a^A{<r)ahf(z)ci^i  +  (o^i)*?*  =  0  (50) 

or 

Now  dehne  quantity  A  such  that 


A  —  a2ia\if  (z)  A{(r)  —  1  (52) 

By  substitution  of  Equation  (52)  into  (51)  and  combining  the  result  with  inequality  (49),  the 
inequality  can  be  written  as: 


2A-t-l 
A  +  1 


>  1 


(53) 
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Thus  the  following  situations  arise: 


j 

A  >  0  =>  dilatancy 
A  <  0  contraction 

A  =  0  =>  constant-volume  deformation.  '  m 


It  is  noteworthy  that  in  deriving  the  above  result,  Equation  (51)  was  written  in  the  form  of 
EJquation  (54). 


2A  = 


(54) 


Thus  the  sign  of  A  is  aiso  determined  by  Ekjuation  (54).  In  other  words,  in  q  -  c  space  there 
exists  a  strught  line  given  by  the  relation 


a^iq  =  a  (55) 

such  that  stress  points  (always  on  or  below  the  failure  envelope)  above  the  straight  line  give  rise 
to  dilatancy  but  stress  points  below  the  line  give  rise  to  contraction.  Points  on  the  strught  line 
give  rise  to  deformation  at  constant  volume.  The  situation  is  illustrated  in  Figure  4. 

Calculation  of  the  Volumetric  Response 

20.  In  this  section,  the  volumetric  response  during  shear  and  its  consequences  are  con¬ 
sidered  in  more  detail  beginning  with  an  explicit  calculation  of  the  intrinsic  time  by  making  use 
of  Elquations  (46)  and  (47).  Substitution  of  the  former  in  the  latter  gives  the  relation 

I/faP  II*  (feP 

■*1 

where  use  was  made  of  Elquation  (23).  Use  of  Elquations  (39)  and  (52)  gives  the  additional  relation 

oiioji  =  0^(1  +  A)  (57) 

At  this  point,  Equations  (56)  and  (5)  combine  to  give  the  following  quadratic  relation  for  de’*/dz. 
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which  gives  two  solutions  for  kde^  fdz.  The  one  is  given  by  Equation  (59) 

=  (59) 

1  +  A  ^  ’ 

The  second  is  equal  to  unity  which  must  be  rejected  since  it  gives  rise  to  the  consequence  that 


de'’ 


dz 


=  0 


for  all  histories  (in  view  of  Equation  (56)),  with  the  implication  that  the  shear  response  is  always 
elastic.  This,  of  course  is  not  admissible,  except  in  situations  where  the  deformation  is  purely 
hydrostatic.  Equation  (59)  gives  rise  to  the  conclusion  arrived  at  earlier  that  if  A  >  0,  then 
de^  <  0  and  the  deformation  is  dilatant,  while  if  A  <  0  the  material  response  will  be  contractive. 
It  follows  from  Equation  (56) ,  in  the  light  of  Equation  (59)  that 


deP 


dz 


1  +  2A 
(1  +  A)2 


(60) 


thus 


dz  = 


1  + A 


\\d^\\ 


(61) 


y/l  +  2A 

Since  dz  is  real  we  require  that  the  radical  should  be  positive,  whereby  the  following  inequality 
must  apply: 


14 


A>-§  (62) 

Recalling  Equationa  (39)  and  (52)  the  implicati<Hi  of  inequality  (62)  is  that,  for  <r  >  0  (compres¬ 
sion) 


2oiiOii  >  (®3) 

Of  particular  importance  is  the  situation  where  the  equality  (64)  applies: 


2oiia$i  =  a  (64) 

in  this  case  and  in  view  of  the  yield  condition  given  by  Equation  (31),  it  follows  that  q  =  0.  Thus 
de^  is  also  zero  and  since  1  +  2 A  =  0,dz  is  indeterminate  from  Equation  (61).  Thus  the  second 
root  of  the  quadratic  Equation  (58)  must  be  chosen,  i.e.,  under  purely  compressive  hydrostatic 
conditions 


as  it  should  be  in  view  of  Bkiuation  (5).  Thus,  the  theory  is  internally  consistent. 


(65) 


Coupling  Mechanisms  Re-examined 

21.  This  section  deals  in  greater  detail  with  the  coupling  mechanisms  that  underlie  the 
deviatoric  and  hydrostatic  stress  responses.  From  a  thermodynamic  viewpoint  the  coefficient  022 
in  the  rate  Ekiuati<m  (2)  accounts  for  the  frictional  dissipation  associated  with  volumetric  plastic 
deformation.  The  indeterminate  hydrostatic  response,  given  by  the  simple  model  under  purely 
hydrostatic  conditions,  is  a  consequence  of  E^quation  (19)  in  which  022  is  proportional  to  a.  The 
model  thus  describes  an  ideal  locking  material  which  displays  no  plastic  volumetric  deformation 
in  response  to  an  increase  in  hydrostatic  stress.  Although  Elquation  (23)  permits  the  introduction 
of  an  independent  hydrostatic  constitutive  equation  there  is  no  apparent  physical  motivation  to 
guide  the  choice  of  such  an  equation. 

22.  In  reality,  the  coefficient  022  in  a  frictional  materied  will  depend  very  strongly  on 
the  prevailing  hydrostatic  stress  because  the  interparticle  frictional  forces  depend  on  the  normal 
forces,  which  are  represented  in  an  average  sense  by  a.  In  addition,  however,  a  signiHcant  part 
of  the  dissipation  under  hydrostatic  conditions  is  due  to  particle  crushing.  It  is  reasonable  to 
suppose  that  the  degree  of  crushing  will  increase  with  the  number  of  contact  points  associated 
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with  an  increase  in  the  density  of  a  soil.  These  two  mechanisms  are  essentially  independent  since 
the  number  of  friction  points  will  not  affect  the  total  frictional  force  (Amonton’s  Law). 

23.  Therefore,  from  a  physical  standpoint,  022  will  depend  on  a  as  well  as  the  porosity 
(or  equivalently  the  density  or  plastic  volumetric  strain)  and  the  dependence  is  likely  to  be 
substantially  additive.  Thus  set 

022  =  (1  - 

where  4>o  is  a  function  of  the  porosity  n,  and  ^1  is  a  constant  such  at  0  <  <  1 

hydrostatic  conditions 

=  ((1  -  4>i)4>oin)  +  ^i<r)  (67) 

which,  in  view  of  the  definition  of  intrinsic  time,  gives  the  following  hydrostatic  stress-plastic 
volumetric  strain  relation,  under  monotonic  conditions  of  continuous  compression: 

Note  that  in  the  case  where  k  is  equal  to  4it  hydrostatic  stress  becomes  infinite.  This  behavior 
is  that  of  an  ideal  locking  material  discussed  previously.  However,  it  is  essential  that  ib  >  to 
ensure  a  positive  (compressive)  volumetric  strain. 

24.  It  is  useful  at  this  point  to  introduce  a  porosity-dependent  function  ^e(n)  such  that 

<(>0  =  (69) 

The  meaning  of  ^«(n)  becomes  evident  in  a  purely  hydrostatic  monotonically  compressive  test 
whereby  Equations  (68)  and  (69)  give 

or  =  4>t{n)  (70) 

Thus,  4>e{n)  is  the  hydrostatic  response  function  under  purely  hydrostatic  conditions  of  increasing 
compressive  stress.  It  is  mentioned  in  passing  that  Equation  (70)  allows  a  connection  with 
Hvorslev’s  (1937)  equivalent  pressure  concept  wherein  (j>t  is  identical  to  the  limiting  hydrostatic 
stress  introduced  in  Equation  (42).  As  in  Hvorslev’s  theory  the  shear  strength  of  the  soil  becomes 
a  function  of  the  ratio 

25.  By  substitution  of  Equation  (66)  into  Elquation  (2)  and  use  of  Equation  (5)  the  yield 
condition  can  be  written  in  the  form  shown  below. 


(71) 


1  ( _ — _ V  =  a* 

\  <*11  J  \(1  —  +  ^i<*^/ 


9’  +  «in*-— -I  + 

It  can  be  readily  shown  that,  if  =  0  and  k  =  Equation  (32)  is  recovered.  In  view  of  Elquation 
(68),  it  is  possible  to  make  Equation  (71)  approximate  the  ideal  frictional  response  of  Equation 
(S2)  as  closely  as  desired  and  still  retain  the  hydrostatic  response  of  Equation  (68).  SpecificsJly, 
we  make  use  of  Equation  (69)  and  define  the  coupling  parameter  c  as 


c  = 


^1 


to  rewrite  the  yield  relationship  as 


9*  + 


(72) 


(73) 


1)^«  +  > 

Regardless  of  the  magnitude  of  Elquation  (73)  may  be  made  as  close  to  Equation  (32)  as 
desired  by  setting  c  as  close  to  unity  as  necessary.  Thus  Equation  (73)  encompasses  the  behavior 
of  idealized  granular  material  while  eliminating  the  indeterminacy  of  the  response  for  purely 
hydrostatic  conditions. 

26.  It  may  be  shown  by  direct  substitution  that,  for  9  =  0,  E}quation  (73)  reduces  to 
Equation  (68).  On  the  other  hand,  the  complete  yield  equation  which  is  analogous  to  Equation 
(45),  can  be  written  as 


qo  ^*V  \»T*/  ^ 


(74) 


where 


C  =  1  + 


1  +  2B>(1  -  B») 


m»(l/m)W4. 
(m(l/m)"/^*B*  -  1)* 


(75) 


and 


B  = 


CO 


(c  -  l)(f>e  +  <^ 

It  is  readily  verified  that  Equation  (74)  is  identical  to  Ekjuation  (45)  when  c  =  1. 


(76) 


27.  Thus  fw,  the  parameter  c  has  been  treated  as  a  device  to  formally  introduce  a 
hydrostatic  yield  stress  into  the  theory;  in  fact,  c  has  an  important  physical  interpretation.  By 
combining  Equations  (27),  (28),  (39)  and  (72)  and  noting  gj,  =  0  the  relationship  for  the  rate  of 
volume  change  is  found  to  be 
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1 _ ^ _ fiV 

The  effect  of  c  >  1  is  to  reduce  the  volume  change  induced  by  shear-volume  coupling  at  small 
values  of  a.  In  view  of  the  assumed  dependence  of  4>e  on  density,  the  degree  that  volume  change 
is  suppressed  by  c  depends  on  the  compactness  of  the  soil.  There  is  experimental  justification  for 
both  the  reduced  shear-normal  coupling  at  low  stress  and  the  influence  of  density  eis  depicted  in 
Equation  (77)  although  quzmtitative  comparison  with  data  remain  for  future  research. 


Summary  of  Simple  Model 

28.  The  simple  model  presented  in  this  part  is  limited  to  isotropic  hardening.  However, 
several  important  principles  outlined  in  its  development  will  have  application  to  the  more  general 
theory  that  is  presented  in  the  next  part; 

a.  The  coupling  between  volume  change  and  shear  is  derived  from  both  the  definition 
of  endochronic  time  and  the  coupling  terms  in  the  rate  equation. 

b.  The  indeterminate  response  to  purely  hydrostatic  loading,  typical  of  models  for 
ideally  frictional  materials,  is  removed  by  introducing  a  hydrostatic  yield  stress. 
Physically,  the  hydrostatic  yield  stress  is  a  measure  of  energy  dissipated  as  a  result 
of  particle  crushing. 

c.  A  yield  surface  was  derived  from  the  definition  of  endochronic  time  and  the  rate 
equations.  By  making  the  resistance  coefficient  proportional  to  mean  stress  a 
model  for  an  ideal  Coulomb  material  is  derived.  An  elliptical  yield  surface  results 
from  the  unrealistic  assumption  that  the  shear  stress  at  yield  is  independent  of 
normal  stress  (as  a  result  of  Eiquations  (27)  and  (36))  even  though  the  rate  of 
energy  dissipation  is  proportional  to  a  (as  a  result  of  Equation  (28)). 

29.  Finally,  it  is  important  to  point  out  one  further  tie  to  conventional  theory  of  pleisticity 
for  frictional  materials.  Owing  to  the  presence  of  stress  terms  in  the  resistance  coefficients  (A), 
the  plastic  strain  increment  vectors  are  not  normal  to  the  derived  yield  surface.  Thus,  the  theory 
developed  is  similar  to  a  plasticity  theory  using  a  non-Msociated  flow  rule.  To  impose  normality 
within  the  context  of  the  theory,  it  would  be  necessary  to  propose  a  purely  cohesive  mechanism 
of  resistance  whereby  [A]  would  not  be  functions  of  stress,  in  which  case  the  frictional  nature  of 
the  behavior  would  be  lost. 
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PART  III:  GENERALIZATION  OF  THE  SIMPLE  THEORY 


30.  The  question  that  poses  itself  is  how  to  generalize  the  simple  model  presented  in  the 
previous  sections.  Though  other  possibilities  may  exist,  it  is  apparent  that  two  choices  present 
themselves  at  the  present  juncture.  One  is  an  internal  variable  generalization  of  Equations  (1) 
and  (2)  whereby  the  plastic  strains  and  e**  are  replaced  by  r  internal  variables  q^**)  and 
respectively.  The  r’th  set  of  the  corresponding  equations  reads: 


3q(r) 


dz 


dz 


=  0 


3^, 


■r  aji 


dq(') 


dz 


dqM  '  dz 

and  and  are  given  by  the  following  expressions  respectively: 


(78) 

(79) 


=  (80) 

r=l 

(81) 

r=l 

This  approach  will  be  developed  in  detail  beginning  on  page  20. 

31.  The  other  choice  presents  itself  in  the  light  of  the  stipulation  that  ai2  =  0  in 
Equation  (1)  and  the  relation 


aji  =  a^js  (82) 

in  the  light  of  which  Equation  (2)  becomes: 

32.  Equation  (83)  may  now  be  interpreted  as  meaning  that  there  are  two  causes  of 
plastic  volumetric  strain.  One  is  the  hydrostatic  stress  and  the  other  the  rate  of  deviatoric  plastic 
work.  However,  from  a  physical  point  of  view,  the  deviatoric  plastic  work  is  external  to  the 
hydrostatic  process,  as  far  as  isotropic  materials  are  concerned  at  least,  and  therefore  it  qualifies 
as  a  thermodynamic  force  of  the  first  kind  in  the  sense  of  Valanis  (1983).  The  concept  was 
discussed  in  great  detail  by  Valanis  in  that  reference.  The  thermodynamic  formal  structure  of 
Equation  (83)  then  becomes 
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(84) 


<T  +  R  =  Cl22~7~~ 
dz 

where  R  is  a  thermodynamic  (external]  force  of  the  “first  kind”  such  that 


R  =  ■  —  (85) 

The  generalization  of  Equation  (84)  to  the  case  of  N  internal  variables  is  given  by  Equations  (86) 
and  (87): 


a  = 


dfP 


“T  ^22  j  — 


dz 


where 


»  (r) 


(86) 


(87) 


(88) 


33.  In  other  words,  insofar  as  the  shear-induced  volumetric  strain  is  concerned,  the  effect 
of  the  deviatoric  plastic  work  is  to  apportion  itself  among  the  N  hydrostatic  internal  mechanisms, 
through  the  coefficients  a^2i-  This  second  approach  will  be  developed  in  greater  detail  beginning 
on  page  28. 

Generalized  Theory  I 

34.  Equations  (1)  and  (2)  are  a  particular  case  of  the  more  general  set  of  equations 
derived  from  linear  thermodynamics  of  internal  variables.  To  show  this,  begin  with  the  relations: 


8  = 


a  = 


deP 

d'ifk 

dfP 


(89) 


(90) 


where  'I'd  is  a  quadratic  function  of  and  is  a  quadratic  function  of  fP  and  q^'’^  and 

g(’')  are  respectively  the  deviatoric  and  hydrostatic  component  of  an  internal  variable  such 
that 


aC'-)  =  qfO  +  IgC*-) 

where  I  is  the  identity  tensor  and 


(91) 
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=  9j(e»,q('>)  +  9r„(e»,g('))  (92) 

The  decompositions  into  deviatoric  and  hydrostatic  components  inherent  in  Ekjuations  (89)  to 
(92)  apply  only  to  the  case  where  the  material  in  question  is  isotropic  in  its  initial  state  and  in 
the  case  of  Equation  (92)  which  is  also  quadratic  in  its  arguments.  For  details  see  Valanis  and 
Read  (1980). 

35.  The  condition  of  positive  internal  dissipation  is  given  by  inequality  (93) 


.aW-I-  <0 

dqW  ^  '  dg(r)^ 

(93) 

whenever 

114(011  >0, 14(01  >0 

(94) 

or 

|4(0l  >  0,  llq(0l|  >  0 

(95) 

In  the  event  that  q^’’^ 

and  are  independent,  inequality  (93) 

reduces  to  the  two  separate 

relations 

0.4W<o,|14''>||>o 

(96) 

and 

(97) 

In  linear  thermodynamics  of  isotropic  materials  the  deviatoric  and  hydrostatic  evolutions  equa¬ 
tions  are  not  coupled.  This  does  not  preclude  shear-hydrostatic  coupling  which  is  brought  about 
by  the  definition  of  intrinsic  time  as  shown  in  Valanis  and  Read  (1980  and  1986).  When  the 
coupling  is  brought  about  only  through  intrinsic  time  the  material  under  pressure  will  always 
densify  when  shear  is  applied.  In  Part  II,  it  was  shown  that  the  presence  of  the  coupling  coeffi¬ 
cient  a2i  in  Equation  (2)  gave  rise  to  both  dilation  and  contraction,  with  dilation  occurring  when 
the  inequality 

a^jKslI  >  (7  (98) 
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is  satisfied  (see  Equation  (49)).  Evidently  dilatancy  cannot  occur  when  ah  =  0. 

36.  In  the  spirit  of  Part  II,  therefore,  and  in  reference  to  Equations  (1)  and  (2)  express 
the  evolution  equations  for  the  internal  variables  in  the  form: 


aqM+““  dz 
dqi^)  dz 


dz 


=  0 


=  0 


(99) 

(100) 


Notable  by  their  presence  are  the  coupling  terms  and  However  in  what  follows  the 
possibility  of  a^j^  =  0  for  all  r  will  be  considered,  thus  abandoning  the  Onsager  symmetry 
condition.  This  was  discussed  in  Part  II.  For  the  purpose  of  integrating  Equations  (99)  and  (100) 
we  express  'id  and  in  the  usual  quadratic  form: 


r=l 


(101) 

(102) 


and  following  Part  II,  write  the  coefficients  011,022  ,  and  821  in  the  form: 


(r)  wp  o(r) 


(103) 


4'.’  =  n4'.'' 


(104) 


and 
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o(r) 
Oji  8 


(105) 


where  Fd  is  the  deviatoric  hardening  function  (in  shear)  and  depends  on  the  hydrostatic  stress, 
Ff^  is  the  hydrostatic  hardening  function  which  depends  on  the  hydrostatic  stress  and  density 
(or  porosity  n).  Both  Fd  and  F/,  are  also  likely  to  depend  on  z.  Equations  (99)  and  (100)  in 
conjunction  with  Equations  (89)  and  (90)  and  in  light  of  Equations  (101)  and  (102)  give  after 
standard  analysis  the  following  constitutive  equations 
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(106) 


where 


and 


" = r  ’  ‘’’S''*' + |j/o’*  ■ 

(107) 

.  dz 

(108) 

<!.»=  — 

*  tn 

(109) 

The  kernel  functions  <f>d  ,  <ftk,  'uid  are  given  in  the  usual  fashion  by  sums  of  positive  decaying 
exponential  terms  according  to  Elquations  (110),  (111)  and  (112) 


N 


r=l 

N 


where 


f=l 

r,  = 


Ar  .  ^ 

-  o(r)  -  -  ^o(r)  «  -  ^r^o[r)^o{r) 


(110) 

(111) 

(112) 

(113) 
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The  second  integral  of  Equation  (107)  arises  from  the  coupling  terms  in  Equation  (100). 

37.  In  deriving  Equations  (106)  and  (107)  it  was  assumed  that  all  internal  forces  were 
zero,  that  is 


dql’’) 


»=o 


=  0 


(114) 


dq('i 


=  0 


(115) 


*=o 


for  all  r.  While  it  is  necessary  that  8  =  0  and  <r  =  0  at  ^  =  0,  Equations  (114)  and  (115)  need 
not  be  satisfied  initially.  In  fact,  it  is  readily  shown  from  Elquations  (89),  (90),  (101)  and  (102) 


that 
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N 


N 


=  E- 


d^H 


(116) 


(117) 


r=l 

Therefore  while  the  sum  of  the  internal  forces  must  be  zero  in  the  initial  condition,  they  may  be 
individually  non-zero.  This  is  a  fact  of  some  importance  as  it  allows  for  the  possibility  of  an  initial 
apparent  shear  stress  history  for  a  sample  subjected  to  a  purely  hydrostatic  stress  history.  Such 
an  apparent  sheu  history  could  result,  for  example,  from  hydrostatically  stressing  an  assemblage 
of  puticles  that  have  preferred  contact  orientations  created  during  deposition.  This  possibility 
is  consistent  with  the  experimental  observation  by  Arthur  and  Menzies  (1972)  that  the  mode  of 
deposition  influences  subsequent  response  to  loading.  The  possibility  of  treating  depositionally- 
induced  anisotropy  by  a  non-zero  initial  state, however,  will  not  be  pursued  hence  but  remains  a 
topic  for  future  research.  The  same  remarks  apply  to  the  development  of  Theory  II  in  the  section 
that  follows. 

38.  Equation  (107)  is  too  complex  for  practical  applications.  To  simplify  the  second 
integral  on  the  right  hand  side  of  this  equation,  first  consider  what  conditions  must  apply  for  the 
relation 


to  hold,  where  kr  are  constants.  An  analysis  using  Equations  (99)  and  (101)  show  that  this  is 
possible  only  if 


In  this  event 


(119) 


<t>d  =  (120) 

that  is,  the  heredity  function  in  shear  is  given  by  a  single  exponential  term.  While  this  is  not 
always,  if  ever,  true,  it  is  a  reasonable  approximation  (to  <f>j)  which  allows  a  significant  reduction  in 
the  complexity  of  the  right  hand  side  of  Equation  (107).  Thus  assuming  (119)  to  be  approximately 
true  then  it  follows  from  Equations  (99)  and  (101)  that 


8  =  - 


E 


r=l 


N  fz 

7^1  ^0 


dz' 


(121) 


24 


But 


dqi') 


(122) 


Thus,  in  view  of  Equations  (118),  (121)  and  (122) 

t,  =  (123) 

r=l 

Now  use  of  Equation  (118)  in  Equation  (107)  gives  the  following  simplified  constitutive  equation 
for  the  hydrostatic  response: 

r(zk  -  (124) 

where 


r(«fc)  =  EfcrC.e-'»'‘‘  (125) 

r=l 

Thus  to  summarize,  the  constitutive  description  of  potentially  dilatant  behavior  in  soils  is  given 
by  Equations  (106)  and  (124)  in  terms  of  two  heredity  functions  <f>d  and  4>h>  the  coupling  ker¬ 
nel  function  r(zh)  and  the  hardening  functions  Fj  and  Fh  in  addition  to  the  elastic  shear  and 
bulk  stiffness  constants  fi  and  K.  The  remainder  of  this  section  considers  the  conditions  where 
Equation  (124)  does  indeed  give  rise  to  dilatant  behavior  in  soils. 

39.  The  dilatant  behavior  of  Equation  (124)  is  studied  most  easily  (and  in  closed  form) 
in  the  particular  case  where  both  4>h  and  F  are  given  in  terms  of  a  single  exponential  function, 
that  is,  when 


(126) 


and 


r  = 


(127) 
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The  multiplier  P  on  the  right  hand  side  of  Elquations  (126)  and  (127)  is  not  put  there  artificially 
but  it  serves  a  useful  purpose  in  that,  in  the  limiting  case  of  very  high  P  ,  both  <f>h  and  F 
tend  asymptotically  to  Dirac  delta  functions,  which  are  useful  approximations  to  and  F  for 
non-cohesive  materials  such  as  sands. 

40.  Using  Laplace  IVansforms  one  may,  in  the  light  of  Equations  (126)  and  (127),  solve 
Equation  (124)  in  terms  of  to  obtain  the  following  relation: 


(128) 


Now  let  <T  be  equal  to  (Tg,  the  constant  hydrostatic  stress  at  which  shearing  begins,  and  let 
dcj  be  the  increment  in  volumetric  strain  due  to  shear  at  constant  hydrostatic  stress.  Then, 
differentiating  Equation  (128)  we  obtain  the  following  relation  for  d^: 


-  F<,K^  dZH 

where  the  following  relation  was  used: 


(129) 


e”  =  £?  +  (130) 

cj  being  the  plastic  volumetric  strain  induced  by  the  initial  hydrostatic  stress  ag.  Thus  in  this 
more  general  approach  the  criterion  for  dilatant  behavior  is  given  by  inequality  (131): 

llsf  >  (131) 

A  0 

41.  The  case  of  a  cylindrical  triaxial  test  at  constant  lateral  pressure  is  more  complex 
since  a  does  not  remain  constant  during  the  test.  Thus  upon  differentiation  of  Equation  (128) 
one  does  not  obtain  Equation  (129). 

42.  Given  that  for  cohesionless  soils  is  substantially  proportional  to  a  ,  then  in  this 
particular  case  and  for  a  >  0  (tr  positive  in  compression) : 


Fd  =  Fla  (132) 

In  view  of  Equation  (132)  the  criterion  for  dilatancy  in  non-cohesive  soils  is  given  by  inequality 
(133) 


i 


I 


I 
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Note  that  condition  (131)  and  its  corollary  (133)  are  identical  in  form  to  condition  (49)  and  (54) 
which  was  obtained  in  the  case  of  the  simpler  theory  of  Part  II. 

43.  In  shear  experiments  on  soils  under  constant  hydrostatic  stress  what  is  actually  mea¬ 
sured  is  the  resulting  volumetric  (plastic)  strain  as  a  function  of  the  shear  strain  (or  shear  stress). 
Thus  what  one  obtains  eventually  is  a  relation  between  the  shear  induced  plastic  volumetric  strain 
and  the  applied  plastic  shear  strain.  To  obtain  such  a  theoretical  relation  from  Ek]uation  (129) 
use  liquations  (5)  and  (109)  to  obtain  the  following  result  on  the  basis  that  dt^  >  0: 


deP 


Oa  Po 


(134) 


(135) 


For  soils  the  hardening  function  F\  is  adequately  represented  by  the  form 

Fa  =  (l-^jl)e^'  +  §a 
*Ph 

where  and  <f>\  are  material  properties  and  has  the  same  definition  as  when  introduced  before. 

44.  Equation  (134)  may  be  integrated  numerically  to  determine  as  a  function  of  e^ 

provided  that  the  material  constants  and  ^  are  known.  Actually  five  material 

constants  are  necessary  to  determine  theoretically  the  shear-induced  volumetric  response  since 
Equations  (132)  and  (134)  dictate  that  and  appear  in  the  ratio  To/Fi. 

Further  Simplification  of  Theory  I 

45.  The  hydrostatic  stress-strain  curve  of  soil  (obtained  under  monotonic  conditions 

of  increasing  hydrostatic  stress)  consists  of  an  initially  convex  part  which  quickly  gives  way  to  a 
concave  part  which  prevails  over  the  remainder  of  the  curve.  The  concave  part  of  the  stress-strain 
curve  is  obtained  readily  by  setting  <t>h{zh)  proportional  to  a  Dirac  delta  function,  i.e.  ,  by  setting 
in  Equations  (126)  and  (127)  =  oo.  In  this  event 


<t>h  =  <f>hS{zk) 


(136) 


r  =  r„«(zfc)  (137) 

46.  With  reference  to  the  purely  hydrostatic  response  (s  =  0)  and  in  view  of  Equations 
(5),  (107),  (124),  and  (135),  it  follows  that; 


s 


\ 


-4 


4 
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r 


1 


a  =  <f>%e^'  (138) 

an  expression  which  agrees  remarkably  well  with  experiment  in  the  concave  region  of  the  hydro¬ 
static  stress-stredn  curve.  Thus 


(To  = 


P 

o 


Thus  in  view  of  Equations  (132),  (135)  and  (139)  Elquation  (134)  becomes: 


(139) 


where 


(140) 


(141) 


This  equation  may  now  be  integrated  numerically  provided  that  the  constants  k  ,  (f>\  and 

To/ FI  are  known. 


Generalized  Theory  II 

47.  In  this  section  the  position  is  considered  that  whereas  the  deviatoric  plastic  work  is 
a  cause  of  hydrostatic  plastic  strain,  it  is  external  to  the  hydrostatic  process  in  the  sense  that 
it  is  not  a  hydrostatic  mechanism.  Therefore,  it  qualiBes  as  a  thermodynamic  internal  force  of 
the  first  kind  in  the  sense  of  Valanis  (1983).  In  the  presence  of  such  forces  the  thermodynamic 
equations  appropriate  to  the  rigid  pleatic  solid  that  represents  the  plastic  behavior  of  the  soil  are 
the  following: 


N 


dqi’") 


(142) 

r=l 

N 

(143) 

r=l 

d^d 
®  deP 

(144) 

°  deP 

(145) 

(146) 
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(147) 


where 


dqM 


dz 


=  Rr 


Rr  =  -4'}a~  (148) 

In  addition,  the  resistance  coefficients  and  022^  are  not  constant  but  are  related  to  the  hard¬ 
ening  functions  Fj  and  respectively  by  the  relations 


where 


Fja 


o(r) 

11  » 


ag  = 


(149) 


Fi  =  Fi{a,z),  F^  =  F^(<T,e’>)  (150) 

48.  The  dependence  of  F4  on  z  is  weak  so  that  for  deformations  other  than  cyclic  histories 
z  plays  a  minor  role  in  Fj  and  may  be  ignored.  Also  the  dependence  of  F/,  on  <t  and  e’’  is  additive, 
a  fact  that  was  discussed  at  some  length  on  page  15.  A  straightforward  analysis  using  Equations 
(142)  to  (148)  in  the  light  of  the  initial  conditions 


q('^(0)  =  0,  ?<’-)(0)  =  0 

gives  rise  to  the  following  set  of  two  constitutive  equations; 


(151) 


where 


(152) 

r  -  *')  + r  ■  *'>’  • 

(153) 

(154) 

,  dz 

(155) 

(156) 

r=l 
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49.  The  determination  of  the  form  of  the  functions  ^d,4>h  &nd  F  from  appropriate 
experiments  is  a  subject  for  future  research.  Some  simplifications  do  occur  when  the  effect  of  the 
deviatoric  plastic  work  rate  is  distributed  uniformly  among  the  hydrostatic  mechanisms,  in  which 
case 

=  all  (for  all  r)  (161) 

50.  In  this  event,  and  in  view  of  Equations  (157)  and  (158) 

r(^fc)  =  a\Mzh)  (162) 

so  that  once  4^kizh)  is  known,  from  a  simple  monotonic  hydrostatic  test,  then  r(2;,)  is  also 
determined  to  within  a  multiplicative  constant. 

Simple  Example  of  Theory  II 

51.  As  an  illustration,  the  theory  is  applied  to  the  simple  case  of  a  cylindrical  triaxial 

test  when  the  soil  is  sheared  under  constant  hydrostatic  stress.  Specifically  the  soil  is  stressed 
initially  under  purely  hydrostatic  conditions  until  a  is  equal  to  (To-  With  a  held  at  this  value  the 
soil  is  then  compressed  axially  and  the  shear  stress,  the  shear  strain  and  the  volumetric  strain 
are  measured.  Also  it  is  assumed  that  <To  is  sufficiently  large  to  lie  on  the  concave  part  of  the 

hydrostatic  stress-strain  curve,  so  that  <t>h{zh)  can  then  be  represented  by  delta  function.  See 

Valanis  and  Read  (1980). 

52.  The  mathematical  consequences  of  this  representation  are  of  interest  in  view  of 
Equations  (157)  and  (160) 


=  (163) 

r=l 

The  limiting  process  by  which  4>h  tends  to  a  delta  function  is  the  one  whereby  fir  -*  which 

event 


=  =  (164) 

r=l 

Of  importance  is  the  fact  that  when  tends  to  a  delta  function,  so  does  r(zji)  by  virtue  of 
Equation  (158).  In  effect 


/J, 


lim  T{zh)  =  ^4^«(^k) 

>r— OO  Fit  “ 


In  this  event 


where 


r(zfc)  =  TcSiznyP^ 


r=l 


The  constitutive  Equation  (153)  now  reduces  the  simple  form 


^  df"  „  deP 

<T  =  (l>g- - 1-  FoB  •  -J- 

dzh  dz 


(165) 


(166) 


(167) 


(168) 


where  Equation  (155)  was  also  used.  This  last  equation  is  used  to  investigate  the  dilateuit  behavior 
of  soils  under  constant  hydrostatic  stress. 

53.  Following  the  deliberations  on  page  18  we  set 


rn 


(169) 


so  that  Equation  (168)  becomes 


dfP  //oP 

<T  =  fc((l  -  <f>H)<t>he^"'  +  ■  d7 

An  analysis  of  the  cylindrical  triaxial  test  shows  that  if  one  sets 


(170) 


»  =  -  ^^s) 


(171) 
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(172) 


where  “1"  denotes  the  axial  direction  and  “3”  one  of  the  lateral  directions,  then 


deP  •  deP  =  (deP)*  (173) 

s  •  deP  =  fldeP  (174) 

Thus  the  equation  that  defines  the  intrinsic  time  z,  on  one  hand,  and  Equation  (170)  on  the  other 
become 

dz*  =  (dcP)*  +  fc*(dfP)*  (175) 


a  =  k((l  - 


(176) 


The  ultimate  goal  is  to  solve  Equations  (175)  and  (176)  simultaneously  and  thus  obtain  a  relation 
between  fP  and  eP. 

54.  To  this  end  note  that  under  purely  hydrostatic  conditions  (s  =  0,eP  =  0)  for 
monotonic  loeuling 


<ro  =  (177) 

for  all  <To.  Therefore  define,  as  before,  a  deviatorically  induced  hydrostatic  strain  such  that 


c5  =  cP  -  eP  (178) 

Equations  (175)  and  (176)  may  now  be  solved  simultaneously  under  conditions  of  monotonically 
increasing  cP  to  give  the  following  result:^ 


der  hr, a-  +  V*'  +  T’s'*  -  1 

where 


(179) 


h  =  (l-^]‘)e^5  +  ^i  (180) 

'Equation  (179)  corresponds  to  the  case  where  s  is  increasing.  A  second  solution  exists  for  Equations  (175)  and 
(176)  which  corresponds  to  unloading.  The  solution  procedure  leading  to  Equation  (179)  is  discussed  in  detail  in 
Appendix  A. 


and 


(181) 


Comparison  of  Theories 

55.  The  theories  presented  in  the  two  previous  sections  differ  as  a  result  of  the  interpre¬ 
tation  of  the  coupling  terms  in  Equations  (78)  and  (79).  Theory  I  results  from  eissuming  that 
coupling  occurs  for  each  mechanism  whereas  Theory  II  results  from  assuming  that  the  coupling 
occurs  through  some  averaging  mechanism  represented  by  the  external  force  R.  Theory  II  is  quite 
amenable  to  simplification  without  loss  of  physical  relevance  whereas  Theory  I  leads  to  relation¬ 
ships  that  are  virtually  intractable  except  for  a  limited  physically  unrealistic  case  in  which  the 
shear  response  is  modeled  by  one  internal  variable.  It  is  worth  considering,  therefore,  how  the 
theories  differ  on  physical  grounds. 

56.  Both  theories  predict  a  contractive-dilatant  response  for  monotonic  loading,  while 
for  unloading  some  important  differences  arise.  Consider  an  experiment  whereby  a  cylindrical 
specimen  is  first  loaded  axially  until  dilation  is  observed  then  unloaded  along  the  same  path. 
In  accordance  with  Equation  (129)  Theory  I  predicts  that  the  volumetric  strain  upon  unloading 
would  be  dilative  arid  would  remain  so  provided  inequality  (131)  was  satisfied.  This  result  is,  in 
fact,  easily  extended  to  the  general  case  described  by  Equation  (107)  by  noting  that  the  coupling 
described  by  the  second  integral  is  a  function  of  the  internal  state  and  independent  of  the  loading 
direction. 

57.  By  contrast  Theory  II  predicts  that  the  sign  of  the  volumetric  strain  increment 
caused  by  coupling  is  controlled  by  the  product  s  •  de**  as  seen  from  either  Equation  (153)  or 
(170).  Thus,  Theory  II  predicts  that  upon  unloading  the  volumetric  rate  should  change  from 
dilation  to  contraction  (see  Appendix  A).  One  of  the  most  important  features  of  the  response  of 
frictional  materials  to  cyclic  loading  is  their  tendency  to  densify  upon  load  reversal.  Therefore, 
Theory  I  appears  to  be  physically  inadmissible  on  this  particular  basis. 


Yield  Surface 

58.  A  feature  of  the  simple  theory  presented  in  Part  II  was  the  existence  of  a  yield  surface 
which  provided  a  means  for  comparison  with  models  developed  from  the  theory  of  plasticity.  For 
multiple  internal  variable  formulations,  such  as  those  presented  in  this  Part,  the  boundary  in 
stress  space  between  elastic  and  plastic  response  is  not  necessarily  marker  by  a  yield  surface.  As 
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shown  by  Valanis  (1980)  for  metals,  a  yield  surface  is  a  mathematical  consequence  of  representing 
at  least  one  of  the  exponential  terms  of  the  memory  kernel  with  a  delta  function.  The  situation 
for  the  models  presented  here  is  similar  as  will  be  illustrated  for  Theory  II.  To  begin,  the  first 
terms  in  the  series  given  in  Equations  (156)  to  (158)  are  approximated  by  delta  functions  or 
alternatively  by  letting  and  0i  become  infinite  while  keeping  the  ratios  Aijai  and 

B\l finite.  In  this  case  the  integrals  in  Equations  (152)  and  (153)  become  Equations  (182)  and 
(183): 

dpP  r**  deP 

8  =  -  z')^dz'  (182) 


dfP  deP 


dfP  deP 


(183) 


which  give  rise  to  a  theory  with  a  yield  surface. 

59.  Equations  (182)  and  (183)  may  be  used  to  obtain  explicit  relationships  for  the  plastic 
strain  rates  in  terms  of  the  stresses  and  history  integrals.  Relationships  analogous  to  Equations 
(46)  and  (47)  (with  ajj  =  0)  can  be  now  derived  £is: 


and 


where 


and 


deP 

q.  = 

(184) 

,<,dfP  deP 

(185) 

(186) 

dfP  /■**  deP 

(187) 

60.  Through  an  analysis  similar  to  that  on  page  15  whereby  the  pleistic  strain  rates  are 
substituted  into  the  expression  for  endochronic  time  the  following  yield  condition  is  obtained: 


8-(s-Q,)  -4>fF^  =  0 


(188) 


61.  Equation  (188)  describes  a  yield  behavior  that  includes  both  kinematic  and  isotropic 
hardening.  If  Fd  and  Fh  are  constant  the  behavior  is  purely  kinematic  hardening.  In  the  case 
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of  an  incompressible  material  where  <l>%Fh  — ►  oo  the  deviator ic  projection  is  circular  as  for  the 
case  of  metals  derived  by  Valanis(l980).  For  compressible  materials  with  Fq  ^  0  the  projection 
is  non-circular  as  a  result  of  the  product  s  •  (s  —  Q«). 

62.  This  is  a  historical  development  of  some  relevance.  Note  that  the  constitutive  relations 
(182)  and  (183)  in  conjunction  with  the  yield  equation  (188)  contain  a  kinematic-cum-isotropic 
hardening  rule  in  the  co-joint  hydrostatic-deviatoric  stress  space.  This  comes  about  naturally 
through  the  vehicle  of  thermodynamics  which  has  not  been  available  to  treatments  that  begin 
with  the  concept  of  the  yield  surfaces. 


PART  IV:  EXAMPLE  APPLICATION  TO  SAND 


63.  In  this  psirt  the  theory  described  on  page  28  will  be  used  to  model  the  response  of 
Sacramento  River  sand  in  the  triaxial  compression  test.  This  test  is  performed  by  increasing  the 
axial  stress  <ri  while  keeping  the  lateral  stress  as  constant.  Thus  the  average  stress  a  increases 
during  the  test  which  complicates  the  analysis  of  the  constitutive  equations.  Therefore,  an  im¬ 
portant  goal  of  the  analysis  in  this  Part  is  to  demonstrate  how  the  theory  can  be  applied  to  the 
more  common  triaxial  test,  an  essential  step  in  use  of  the  model  for  practical  applications.  The 
data  used  for  the  analysis  is  taken  from  Lee  and  Seed  (1967).  These  data  are  presented  in  Figures 
5  and  6. 

Hydrostatic  Response 

64.  The  hydrostatic  response  is  described  by  the  one-term  approximation  presented  by 
Equation  (168).  To  simplify  the  determination  of  parameters  a  form  of  is  used  to  yield  a 
relationship  equivalent  to  Equation  (169)  given  by^: 

<f>oFh  =  -  4>i)<l>e{n)  +  <f>ia  (189) 

Note  that  <f>i  and  ^e(n)  were  previously  defined  with  regard  to  the  simple  model  described  in 
Part  II  (see  page  15)  where  0  >  >  1  and  <f>e{n)  is  a  density-dependent  hydrostatic  limit  stress. 

Also  note  that  because  <f>o  always  appears  as  a  multiplier  to  F^,  we  can  set  (po  =  1  without  loss 
of  generality.  For  purely  hydrostatic  loading  Equation  (168)  reduces  to: 

a  =  ^*(n)  (190) 

65.  For  consistency  with  the  data  presentation  in  Figure  5  the  specific  volume  will  be 
used  which  is  related  to  the  porosity  by: 


V  =  —  (191) 

Also,  for  purposes  of  computation  it  is  desirable  to  express  ao{v)  in  terms  of  plastic  hydrostatic 
strain.  The  change  in  hydrostatic  strain  can  be  related  to  the  change  in  void  ratio  using  the  well 
known  relationship: 

‘Note  that  these  relationships  differ  from  those  in  Equations  (66)  to  (70)  because  of  the  appearance  of  k  in  the 
definition  of  dzi,. 
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1.90 


Figure  5:  Pressure-specific  volume  curves  for  Sacramento  River  sand  at  four  initial  densities 


Oi/Oa 


Figure  6:  Stress-strain- volume  change  data  for  Sacramento  River  sand  with  Vo  —  1.87  and 
<T3  =  2kg/cm* 
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It  will  be  assumed  that  the  initial  unstressed  state  of  the  sand  corresponds  to  the  loosest  state. 
The  limiting  pressure,  can  be  expressed  as  a  function  of  specific  volume  as 

if>e  =  (193) 

where  Vg  is  the  loosest  state  obtainable  by  the  soil.  This  state  is  taken  as  the  so-called  maximum 
void  ratio  (cmaz)  which  defines  the  loosest  packing  that  can  be  achieved  by  a  standardized  test. 
Lee  and  Seed  reported  e^a*  for  Sacramento  River  sand  to  be  1.003.  Thus  Vg  =  2.003. 

66.  Equation  (193)  implies  that  the  initial  state  does  not  depend  on  details  of  hydrostatic 
loading  history.  For  example,  a  loose  sand  and  dense  sand,  if  loaded  to  sufficiently  high  stress, 
can  be  brought  to  the  same  u,  even  though  the  two  specimens  have  significantly  different  stress 
histories;  as  far  as  mechanical  response  to  subsequent  loading  is  concerned,  the  two  samples  would 
be  identical.  For  purely  hydrostatic  histories  this  simplistic  interpretation  of  the  sample  behavior 
can  be  justified  based  on  the  use  of  Equation  (166)  and  indeed  the  assumed  equivalence  between 
plastic  strain  and  void  ratio  is  implied  in  many  soil  mechanics  models  based  on  normalized 
behavior  (e.  g.  Ladd  and  Foott,  1974).  In  practice  this  approach  may  introduce  significant  error 
because  it  ignores  the  apparent  shear  history  introduced  through  sample  densification.  The  role 
of  sample  preparation  in  the  context  of  the  initial  condition  was  discussed  in  Part  III,  page  23. 
Nevertheless,  it  will  be  assumed  here  that  all  specimens  are  initially  at  a  state  that  can  be 
described  by  an  equivalent  purely  hydrostatic  loading  history  that  is  specified  by  the  scalar  v. 


Shear- Volume  Coupling 

67.  Evaluation  of  the  coupling  between  shear  and  volume  change  begins  with  Equa¬ 
tion  (168)  which  can  be  written  for  the  triaxial  compression  test  conditions  through  substitution 
of  Equations  (109)  and  (171)  as: 

deP  deP 

,  =  (194) 

By  combination  of  Equations  (175)  into  (194)  the  following  relationship  for  rate  of  volume  change 
is  derived  in  Appendix  A: 


deP 


1-r*— 


T, 

r<,-V  + 

a* 


hi 


(195) 


Equation  (195)  is  valid  for  increasing  s  and  in  contrast  to  Equation  (179)  is  not  restricted  to 
shear  histories  with  constant  hydrostatic  stress.  Note  that  the  condition  for  dilatancy  is  identical 
to  that  derived  on  page  28  for  constant  a. 

68.  The  two  constants  To  and  k  can  be  determined  from  the  slope  of  the  a  plot  of  e’’ 
versus  e^.  The  constant  is  the  ratio  s/<7  at  the  point  in  the  test  where  de^  =  0  (initiation 
of  dilation).  The  coupling  constant  k  is  computed  by  substitution  of  measured  dilatancy  rates 
dt^ jdf?’  into  Equation  (195).  Prom  Figure  6,  d^ JdtP  =  0  at  aja  =  1.05  which  gives  F,,  =  0.95. 
At  aja  =  \.\A,d€^ ! del'  —  -0.056  and  Fhja  =  2.97.  By  tried  computation  4>x  =  .8  wea  found  to 
give  a  good  fit  to  the  data.  Substitution  of  these  values  into  Equation  (195)  gives  k  =  0.492. 


Shear  Response 

69.  The  shear  response  for  the  triaxial  compression  test,  obtained  by  substitution  of 
Equations  (171)  and  (172)  into  Equation  (152),  is  given  as  follows: 


(196) 

The  above  equation  can  be  written  as: 

3  =  f  4>d{zd-  z')G[z^)dz' 

Jo 

(197) 

where 

(198) 

Owing  to  the  small  values  of  k  and  noting  de^fd^  <  1  it  is  readily  determined  from  Equation 
(175)  that  de^/dz  «  1.  Therefore,  evaluation  of  G(zd)  amounts  to  determining  Fd. 

70.  Equation  (198)  is  complicated  by  the  fact  that  G{zd)  is  proportional  to  a  which, 
in  contrast  to  the  simple  examples  of  Part  III,  varies  over  the  history  of  the  test.  A  simpler 
relationship  can  be  derived  by  noting  that  the  history  of  a  is  tied  to  the  history  of  s  through  the 
following  linear  relationship: 


<7  =  (To  +  5s 


(199) 
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where  <To  is  the  hydrostatic  stress  before  shearing  and  6  is  a  constant  equal  to  zero  for  constant 
hydrostatic  stress  and  .408. . .  in  the  case  of  constant  lateral  stress.  Upon  substitution  of  Equation 
(199)  into  Equation  (197)  a  linear  integral  equation  in  a  is  obtained  which  can  be  solved  through 
a  straight-forward  application  of  Laplace  transforms  once  the  kernel  function  is  defined.  For 
the  example  considered  here  the  kernel  function  will  be  expressed  as  a  delta  function  and  one 
exponential  term  given  by  Equation  (200): 

M^d)  =  (200) 

which,  upon  substitution  into  Equation  (197)  gives: 


=  ffo  ^^d  +  Ai  e‘‘-‘‘dz^  +  h  [4,0^8  +  e^‘-*‘8dz^ 


Solving  for  a  gives: 


(201) 


I 

where 


and 


a  —  Soo  "I"  (®o  Soo)c  * 


= 


<f>i 

l-bcf>i 


Soo  =  (^^ai  +  Ai)^(l^“^^,) 


a  =  ai 


6i4i 

l-bAo 


(202) 

(203) 

(204) 

(205) 


71.  Based  on  the  peak  stress  value  shown  in  Figure  6,  Soo  =  4.28  kg/cm*.  The  two 
constants  So  —  Soo)  and  a  can  be  obtained  from  the  test  data  as  illustrated  in  Figure  7.  Using 
the  data  from  this  figure  it  is  found  that  o  =  88.5.  These  data  along  with  the  elastic  constants 
are  sufficient  to  compute  the  response  of  the  triaxial  compression  test  with  constant  <73.  The 
parameters  needed  for  analysis  under  general  loading  conditions  can  be  computed  from: 


<f>i  = 


So 


(To  +  bag 
(Tg  -f  bag 
(Tg  +  bS, 
Soo 


-a 


Qtl 


(206) 

(207) 

(208) 
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Figure  7:  Semi-logarithmic  plot  of  data  for  determination  of  constants. 
Comparison  to  Experiment 

72.  A  comparison  between  the  theory  and  test  data  are  shown  in  Figure  8.  The  shear 
strain  e  was  computed  from  the  definition  of  plastic  strain  given  by  Equation  (209); 

eP  =  e  -  s/2m  (209) 

where  2/i  =  1700kg/cm*  and  e**  is  found  by  integrating  dc’’  =  Fddzd  in  view  of  Equations  (202) 
and  (199),  and  Fd  =  to  get: 

^  =  {<ro  +  bSoo)zd  -  ^(«oo  -  »o)(l  -  c~“*')  (210) 

Thus  the  shear  response  can  be  expressed  as  two  parametric  equations  in  Zd-  The  volumetric 
response  is  obtained  by  numerically  integrating  Equations  (189),  (192),  (193)  and  (195)  to  obtain 
c^.  The  elastic  volumetric  strains  during  shear  are  computed  using  an  elastic  bulk  modulus 
K  =  1133  kg/cm’. 

73.  The  comparison  shown  in  Figure  8  is  reasonably  good  for  larger  strains.  At  small 
strains,  the  error  created  by  the  one-term  approximation  is  evident,  especially  in  the  initial  shear- 
induced  volume  change  where  the  plastic  component  of  strain  is  underestimated.  Note  that  by 
use  of  more  exponential  terms  in  Equation  (200)  the  initial  elastic  response  can  be  made  as  small 
as  required  to  fit  the  data. 
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Shear  Stress  —  S,  kg/cm 


PART  V:  CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 


74.  The  work  described  in  this  report  wm  intended  to  provide  a  theoretical  framework 
on  which  to  derive  constitutive  models  for  frictional  materials  based  on  linear  thermodynamic 
relationships.  The  principal  goal  of  this  work  was  to  develop  a  fundamentally  sound  model  which 
correctly  predicted  both  contraction  and  dilation  caused  by  shearing.  Two  results  drawn  from 
this  work  are  as  follows: 

a.  The  contractive-dilative  response  of  soil  is  captured  through  a  straight-forward 
application  of  the  theory  of  endochronic  plasticity.  In  particular,  dilatancy  results 
from  coupling  terms  in  the  rate  equations.  A  significant  finding  is  that  coupling 
among  internal  mechanisms  is  through  an  averaging  process  that  can  be  treated  as 
an  thermodynamic  force  which  acts  external  to  the  hydrostatic  mechanism.  This 
is  in  contrast  to  the  equally  plausible  hypothesis  that  there  is  an  independent 
coupling  between  the  shear  and  normal  components  of  each  mechanism. 

b.  The  model  correctly  predicts  that  contraction  occurs  immediately  upon  unloading 
even  for  the  case  where  dilation  occurs  in  the  loading  segment  of  a  stress  path. 
This  is  in  direct  agreement  with  experiments. 

75.  Future  research  should  be  directed  as  follows: 

a.  The  model  should  be  extended  to  the  case  of  undrained  behavior.  Procedurally, 
this  could  be  done  quite  simply  by  applying  the  constraint  de  =  nCu,du  where 

is  the  compressibility  of  pore  fluid  and  du  is  the  increment  of  pore  pressure.  The 
resulting  relationships  could  then  be  used  for  analysis  of  undrained  tests.  Among 
the  goals  of  this  study  would  be 

•  determination  of  properties  from  undrained  tests 

•  prediction  of  liquefaction 

•  comparison  of  soil  response  for  both  drained  and  undrained  conditions. 

b.  A  detailed  study  on  property  determination  is  needed.  In  particular,  the  determi¬ 
nation  of  the  shear-kernel  function  from  data  for  complicated  stress  paths  and 
generalized  forms  of  Fj  requires  investigation. 

c.  The  dependence  of  elastic  constants  ft  and  K  on  stress  and  plastic  strain  requires 
investigation. 

d.  The  response  of  the  model  to  non-proportional  stress  paths  in  deviatoric  stress 
space  and  for  loading  under  rotating  principal  stress  axes  requires  investigation. 
The  kinematic  hardening  behavior  described  by  Equations  (182)  to  (188)  appears 
to  be  consistent  with  experimental  data,  reported  by  Alawaji  et.  al.  (1986).  That 
data,  documented  in  detail  by  Sture,  Alawi  and  Ko  (1988),  should  be  used  to 
investigate  the  kinematic  hardening  mechanism  further. 
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:« 

e.  The  analysis  needs  to  be  extended  to  consider  the  effects  of  the  third  stress  in¬ 
variant. 

f.  The  model  should  be  implemented  in  two  and  three  dimensional  numerical  analysis  ^ 

codes  for  practical  boundary  value  problems. 

g.  The  problem  of  uniqueness  of  solutions  to  bound  ary- value  problems  formulated 
from  the  model  (already  initiated  by  Valanis  (1989))  should  be  studied  further. 
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APPENDIX  A:  DERIVATION  OF  STRESS-DILATANCY  RELATIONSHIP 


1.  In  this  Appendix,  the  derivation  of  Equation  (195)  from  (175)  and  (194)  is  described 
in  detail.  Relationships  such  as  Equations  (134),  (140),  (179),  and  (195)  which  express  the 
dilataney  rate  de^/de^  to  the  sheau^- normal  stress  ratio  s/a  are  referred  to  as  stress-dilatancy 
relationships  in  soil  mechanics  literature.  The  derivation  of  these  relationships  in  the  present 
analysis  is  b2ised  on  combining  the  rate  equations  with  the  definition  of  endochronic  time  and 
solving  the  resulting  simultaneous  equations  for  df** /de**.  Whereas  dc**/ de*”  appears  sis  a  quadratic 
term  in  the  equations,  two  solutions  are  possible.  Therefore,  the  main  task  of  this  appendix  is  to 
describe  how  the  appropriate  solution  is  chosen. 

2.  To  reduce  the  complexity  of  the  relevant  equations,  define  the  following; 


dx  =  kd^ 

(Al) 

dy  —  de” 

(A2) 

Fh 

a  =  — 

a 

(A3) 

b  =:  r.- 

a 

(A4) 

Write  Equations  (194)  and  (175)  respectively  as 

=  > 

dz  dz 

(A5) 

(A6) 

where  the  solution  is  sought  in  terms  of  dx/dy.  These  equations  may  be  interpreted  geometrically 
as  the  intersection  of  a  unit  circle  (Equation  (A6))  by  a  line  (Equation  (A5))  as  plotted  on  the 
cartesian  axes  dy/dz  and  dx/dz  in  Figure  Al.  The  two  solution  points  are  5i  and  S^.  The 
lines  connecting  the  origin  with  these  points, Lj  and  Lj,  have  slopes  dy/dx  which  in  view  of 
Equations  (Al)  and  (A2)  are  the  two  possible  dilataney  rates. 

3.  The  solutions  to  (A5)  and  (A6)  can  be  determined  upon  substitution  of  Equation  (175) 

to  be 


dx  _  1-6* 

dy  ah  ± 

With  reference  to  Figure  Al,  Si  corresponds  to  the  positive  solution  whereas  corresponds  to 
the  negative  solution,  that  is 


Al 


dx 

1-6^ 

for  Si 

(A8) 

dy  ~ 

ab  +  yja}  —  1 

dx 

1-6* 

for  S2 

(A9) 

dy  ~ 

a6  —  \/a*  -f  5*  —  1 

4.  The  line  given  by  Equation  (Al)  crosses  the  dy/dz  and  dx/dz  axes  respectively  at 
points  Pi  and  P2.  The  ordinate  value  of  Pi  is  1/6  and  is  a  measure  of  the  shear-normal  stress 
ratio;  at  s  =  0,  1/6  is  infinite.  The  demarcation  between  contraction  and  dilation  (e.g.  inequality 
(98))  lies  on  the  circle  at  6  =  1;  thus  for  the  loading  condition  dilatancy  occurs  when  Pi  lies  inside 
the  circle. 

5.  Point  Pj  crosses  the  abscissa  at  1/a,  which  is  a  measure  of  the  ratio  between  the  pre¬ 
vailing  hydrostatic  stress  and  the  hydrostatic  yield  stress.  Under  hydrostatic  loading  conditions, 
pleistic  yielding  begins  on  the  circle  at  a  =  1.  In  general,  a  <  1.  Therefore,  P2  lies  on  or  within 
the  circle.  ^ 

6.  When  Pi  lies  above  the  origin,  the  axial  stress  is  greater  than  the  lateral  (confining) 

stress  (commonly  referred  to  as  a  triaxial  compression  test)  and  ai  >  >  0.  When  Pi 

lies  below  the  origin,  the  axial  stress  is  less  than  the  confining  stress  (triaxial  extension)  and 
cri  =  a2  >  (Ta  >  0.  Points  5i  or  S2  falling  above  the  abscissa  correspond  to  solutions  with 
increasing  e’’  whereas  points  falling  below  the  abscissa  correspond  to  decreasing  values  of  e^.  A 
load  reversal  in  a  triaxial  compression  test  will  be  characterized  by  solutions  dy  <  0  with  1/6  >  0; 
that  is,  a  load  reversal  is  depicted  by  Pi  above  the  origin  and  a  solution  point  falling  on  the  circle 
below  the  abscissa. 

7.  The  solution  for  the  initial  condition  (s  =  0)  in  a  triaxial  test  with  a  =  is  shown  in 
Figure  A2.  In  this  case.  Equation  (A5)  describes  a  vertical  line  which  is  tangent  to  the  circle  at 
(l/a,0)  giving  as  the  only  solution  d^ fdz  =  1  (or  equivalently  de^ [de^  —  00).  Two  solutions  are 
possible  for  the  initial  state  when  a  <  1  as  shown  in  Figure  A3.  The  positive  solution  S\  would  be 
chosen  for  a  triaxial  compression  test  (e^  increasing)  while  the  negative  solution  Sj  corresponds 
to  the  triaxial  extension  test  (e^  decreasing). 

8.  Figures  A4  and  A5  illustrate  the  solution  as  s  is  increeised  in  the  triaxial  compression 
test.  For  1/6  >  1,  both  the  positive  and  negative  solutions  correspond  to  contraction  as  both 

*In  soil  mechanics  terminology,  o  =  1  corresponds  to  the  normally  consolidated  state  while  o  <  1  corresponds 
to  an  over  consolidated  state.  A  dense  sand  under  low  confining  pressure  or  a  highly  over  consolidated  clay  would 
have  ft  1 . 


A2 


dy 


Figure  Al:  Solutions  for  stress-dilatancy  equation. 


points  of  intersection  {Si  and  S^)  fall  to  the  right  of  the  origin.  The  positive  solution  corresponds 
to  loading  while  the  negative  solution  applies  upon  unloading  (when  de^  changes  sign).  As  s  is 
increased  such  that  i/b  <  1,  the  positive  (loading)  solution  falls  to  the  left  of  the  origin  which 
corresponds  to  dilation.  Here  again,  the  negative  (unloading)  case  corresponds  to  contraction. 
The  extension  test  is  similarly  analyzed  leading  to  the  conclusion  that  dilation  is  the  result  of 
loading  beyond  |l/6|  <  1,  with  contraction  resulting  otherwise.  The  various  cases  that  may  occur 
in  a  cyclic  load  history  are  illustrated  by  the  sequence  of  sketches  shown  in  Figure  A6. 

9.  Finally,  note  that  by  writing  the  product  s  •  de’’  a3  (Idea’ll  ||8(|  cos  ip,  where  cos  0  =  1  for 
loading  and  coaip  =  —  1  for  unloading,  that  all  cases  can  be  expressed  in  the  form  of  the  solution 
given  by  Equation  (AlO). 
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Figure  A5:  Solutions  for  dilative  case  in  triaxial  compression  test. 


Figure  A6:  Correct  solution  (shown  as  solid  point)  for  various  points  in  compression-extension 
experiment  of  triaxial  test. 


